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We revisit the work of K. Goeke, et al. [[jj who considered a chiral model for the nucleon 
based on the linear sigma model with scalar-isoscalar scalar-isovector mesons coupled to quarks 
and solved using the coherent-pair approximation. In this way the quantum pion field can be 
treated in a non-perturbative fashion. In this work we review this model and the coherent pair 
approximation correcting several errors in the earlier work. We minimize the expectation value of 
the chiral hamiltonian in the ansatz coherent-pair ground state configuration and solve the resulting 
equations for nucleon quantum numbers. We calculate the canonical set of nucleon observables and 
compare with the Hedgehog model and experiment. Using the corrected equations yield slightly 
different values for nucleon observables but do not correct the large virial deviation in the 7r-nucleon 
coupling. Our results therefore do not significantly alter the conclusions of Goeke, et al.. 
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CO . I. INTRODUCTION 

t> 

■ It is widely believed that quantum chromodynamics (QCD) is the fundamental theory underlying the strong inter- 
action. Regrettably, reliable first principle calculations of hadronic structure and reactions based on QCD are still 
00 ■ some time off. Nevertheless, simpler QCD-motivated phenomenological models have been proposed which preserve 
ON ' the important property of chiral symmetry. Familiar examples include Skyrme-Witten models Q and hybrid chiral- 
soliton models such as those of Kahana and Ripka Birse and Banerjee Birse ||, Fiolhais, Goeke, Griimmer, 
O { and Urbano ||, Goeke, Harvey, Griimmer, and Urbano hereafter referred to as GHGU and others (see Ref. 0). 

Such approaches, in particular Ref. Q, argue that spontaneous symmetry breaking of the QCD Lagrangian gives rise 
to an effective chiral Lagrangian of the Gell-Mann-Levy sigma-model form involving explicit quark, scalar-isoscalar 
meson (sigma, a), and pseudoscalar-isovector meson (pion, 7?) degrees of freedom. The model has been solved in 
• • ■ mean-field using the "Hedgehog" ansatz which assumes a configuration-space-isospin correlation for the pion field, 
7T = rir, and for the quarks. One drawback to this ansatz is that it breaks both rotational and isospin invariancc 
(although "grand spin" , K = I + J, remains conserved) requiring some projection onto physical states at the end. 
Considerable attention has been given to the problem of projection in the calculation of observables ||. 

In spite of this drawback the model is quite successful at predicting baryon properties. Constraining the pion mass 
and decay constant with experimental values, the model contains but two additional free parameters (which can be 
written in terms of the effective quark and sigma masses), yet the model makes quite respectable predictions for a 
host of hadronic properties (mass, magnetic moments, sigma commutator, pion-nucleon coupling constant (g^NN), 
axial vector coupling constant {gji), as well as weak and electromagnetic form factors). 

Another criticism of this approach, however, is the approximate treatment of the pion whose very light mass argues 
for a quantum, as opposed to mean-field, treatment. However, treating light strongly-coupled fields is in general 
difficult as perturbative methods are unreliable. In part to address such issues, the coherent-pair approximation was 
developed by Bolsterli J^] and used by GHGU to treat the pions as true quantum fields non-perturbatively. We find 
this approach promising and in this work revisit the hybrid chiral model of Goeke, et al. correcting several errors. 
Besides the inclusion of quantum pionic degrees of freedom non-perturbatively, another advantage is that one need 
not invoke the so-called Hedgehog ansatz for the pion field. Specifically, the symmetry of the quarks with respect to 
permutation induces the space-isospin correlations in the pion field. Following the ground-breaking work of GHGU, 
we minimize the chiral hamiltonian with respect to the coherent-pair Fock-space ansatz ground-state configuration, 
then calculate nucleon properties, and compare with other chiral models. We find that once minimization with respect 
to the coherence parameter is achieved, the results for calculated observables are not significantly different from those 
calculated in GHGU. In addition we present results for the sigma commutator which measures the degree of chiral 
symmetry breaking. 
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For completeness and ease of comparison the organization of this paper closely follows that of GHGU. We present 
the starting Lagrangian and the major intermediate results, refer the reader to the original paper for details, and note 
where our results differ from GHGU. Where appropriate, we will give the corresponding GHGU equation references. 
We present the background theory in section II, a summary of the coherent-pair approximation results in section III, 
the variational equations in section IV, the derived nucleon properties in section V, the numerical results in section 
VI, and our summary and conclusions in section VII. We find that the corrected equations do not significantly alter 
the results originally reported in GHGU. When compared with the Hedgehog model the coherent-pair approximation 
shows systematically smaller mesonic contributions to the nucleon observables and the energy densities which may be 
related to poor self-consistency in the pionic sector. 



II. CHIRAL QUARK-MESON MODEL 

We begin with the chiral model Lagrangian of Gell-Mann and Levy [jlO| but with explicit quark degrees of freedom. 
(Discussion of how such a form may be argued from QCD are given in Refs. [2 and 5].) After chiral symmetry 
breaking, inducing a pion mass, the Lagrangian can be written (GHGU Eq.[2.1]) as 



C{x) = iipftip + ^(d^ad^a + ■ d^ir) + gip(a + ij 5 f ■ if)ip - U(&, tt), (1) 



with 

W^iJ^^^+^-i/ 2 ) 2 -/^- (2) 

where the carat (") denotes a quantum field, /„• is the pion decay constant, is the pion mass, and v, g and A are 
constants to be determined. In the standard scenario spontaneous symmetry breaking generates masses for the quark 
and sigma fields and the linear sigma term, which breaks the chiral symmetry and generates the small pion mass 
which would be zero otherwise as the Goldstone boson of the theory. The vacuum then has a unique non-vanishing 
scalar field expectation value, 

™ = =► fro = 0, ^ = =► a = f n . (3) 

Then the three undetermined constants in the original Lagrangian can be written in terms of the three effective 
masses: m q = —gcro, m 2 , = \ 2 {S<7 2 — v 2 ), and m 2 = \ 2 ((7q — v 2 ). 

We take the experimental values f v = 93 MeV and = 139.6 MeV, leaving g and m a as the only free parameters 
which must be determined. The additional parameters introduced by the coherent-pair approximation are constrained 
by minimization. 

Introducing the conjugate momenta, one formally converts the Lagrangian density to a Hamiltonian density (GHGU 
Eq. [2.10]) into 

H(r) - \{P a {r) 2 + [W(r)] 2 + P n (r) 2 + [W(r)] 2 } 

+U(a, #) + ft (r) (-iaV)V'(r) (4) 
-gft(r)(/3&(r) + i/3~/ 5 r ■ 7?)V>(r), 

where a and f3 are the usual Dirac matrices. In the above expression ip, a, and tt are quantized field operators with 
the appropriate static angular momentum expansions (GHGU Eqs. [2.11], [2.17], and [2.13]), 



*(r) = [|]4 / jM'i^-!— ]i £ j,(fc)Y;„(n r )[4;; >(k) 

_ * ' Pirn in 



+ (-)' 



i+w - 1 — i" 
u t-m 
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$(r) = (( r I n i mw )dl] m + ( r I n 3 mw )dl]J), (7) 

njmw 

where the \njmw) and \njmw) form a complete set of quark and antiquark spinors with angular momentum quantum 
numbers and spin-isospin quantum numbers j, m, and w, respectively. The notation is slightly altered from GHGU 
in that isospin labels appear as superscripts and spin labels appear as subscripts. The corresponding conjugate 
momentum fields have the expansions (GHGU Eqs. [2.14] and [2.18]), 

P„(r) =i^°°d 3 fc[|^] l [c t (fc)e- <fcV -c(fc)e+ ifcV ], (8) 



4(r) =i[|]» J dkk 2 [^]i £ U{kr)YUaZ\k) 

{-) m+w a\Zl{k)]. (9) 

Here the c(fc) destroys a cr-quantum with momentum fc and frequency u> a {k) — (k 2 + m 2 )^ and a\™ m (k) destroys a pion 
with momentum k and corresponding = (k 2 + m 2 )i in isospin-angular momentum state {£m; iw}. For convenience 
one constructs the configuration space pion field functions needed for the subsequent variational treatment by defining 
the alternative basis operators, 

b\Z = j dkk 2 Uk)a\Z{k), (10) 

where £,t{k) is the variational function. Taking this over to configuration space defines the pion field function (GHGU 
Eq. [3.11]), 

In the following only the i = 1 value is used and the angular momentum label will be dropped. 



III. GROUND STATE CONFIGURATION ANSATZ 



The ansatz Fock state for the nucleon is taken to be (GHGU Eq. [4.23]) 

\NT 3 J Z ) = [a(\n) ® |P °))t 3 ./, + P(\n) ® |P 1 1 ))t 3 ./, + 7(1* > ®\P?))t 3 jM), (12) 

where |S) is the coherent sigma field state with the property (GHGU Eq. [3.7]): (E|o-(r)|E) = a(r), and |P °) (|Pf^)) 
are pion coherent-pair states to be determined. The normalization of the nucleon state requires a 2 + (3 2 + j 2 = 1 . 
The permutation symmetric form of the SU(2) x SU(2) x SU(2) quark wavefunctions imply that the source terms 
in the pion field equations will induce an angular momentum-isospin correlation for the pion field. Thus, since the 
pion is isovector, the only allowed angular momentum of the pion will be t — 1, so in the treatment to follow only the 
I = 1 term of the pion field expressions Eqs.iJ and II is retained. 

One constructs a pionic coherent-pair state with quantum numbers of the vacuum as follows Consider the 
scalar-isoscalar coherent state, 

I^HEt^i 1 ^ 1 ™, (13) 

where the double-dot notation refers to spin-isospin (i.e. {m,w}) scalar contractions. A coherent state with spin- 
isospin={ll}, \Pi£l), is constructed by operation of (— ) m+w b\Z m upon \Pq). Sequential such operations appropriately 
contracted are assumed to close, yielding a recurrence relation for f n , namely, 

H m+W b\zZ\P°°) = a\PlZ)- (14) 

A contracted second application gives 



3 



^b\Z\P^) = Lb\P°), 



(15) 



where the spin-isospin multiplicity is L — (2L + 1)(2T + 1) = 9 and where closure is forced by associating the vacuum 
quantum number multipion state uniquely withlPp). The coherence is determined by the parameter, x, defined by 



(b{i:b{i)\P}=x\P). 



(16) 



Here |P) can be either |Pq) or IP 1 ) and x = Lab = 9ab serves as a (as yet free) coherence parameter and the symbol 
b\ t : b\ ^ indicates the coupling to a scalar-isoscalar. Inserting these into Eq.(p3|) and Eq.(16), we obtain the recursion 
relation (GHGU Eq. [4.11]), 



_ x(2n + l) t 



[L + 2n) 



which can be solved to give (GHGU Eq. [4.12]) 



_ a"(2n-l)l!(L-2)ll < 
/n ~ (L-2 + 2n)!! /0 ' 



(17) 



(18) 



where /o is given by the normalization of \Pq). The a and b parameters can be expressed as functions of x from the 
normalization of |Pi). We obtained values for the function / , a(x) and c(x) (b is determined from a) different from 
those of GHGU [1988]. We find 



a(x) 



(105 + 45a; 2 + x 4 ) sinhx - (105 + 10a; 2 )a;cosha; 



— (15 + 6x 2 ) sinha; + (15 + .T 2 )a;cosh: 



(19) 



c(x) 



and 



with y = x 2 and 



-(945 + 420a; 2 + 15a; 4 ) sinhx + (945 + 105a; 2 + a; 4 )a;cosha; 



(105 + 45a; 2 + a; 4 ) sinha; - (105 + 10a; 2 )a;cosha; 



/ - 2 = (i-2)!!2 £ - 1 ar 1 cosh(x), 



(2L + 1) -JL-(-L\JL 
2 ' v ~ dy ~ { 2x'dx' 



(20) 



(21) 



(22) 



which should be contrasted with GHGU Eqs. [4.13], [4.21] and [4.22]. In addition we found that we needed another 
factor for the jP 1 ) matrix element of the four-pion term (implied summation over repeated indices), 



(Pl\blbyb p b a \Pl)=8ld(x) 2 , 
where greek subscripts include both spin and isospin and d(x) is given by 



d(x) = 



(7560 + 3465a; 2 + 165a; 4 + a; 6 ) sinha; - (7560 + 945a; 2 + 18a; 4 )a; cosh a; 



(105 + 45a; 2 + a; 4 ) sinha; - (105 + 10x 2 )x cosh a; 



1/2 



(23) 



(24) 



IV. THE FIELD EQUATIONS 

The total energy of the baryon is given by 

/>oo 

E B =< BT 3 J Z \ / d 3 r : H(r) : \BT 3 J Z >, 



(25) 
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where B = N or A. The field equations are obtained by minimizing the total energy of the baryon with respect 
to variations of the fields, {u(r),w(r) 7 <j(r) 7 <f)(r)}, as well as the Fock-space parameters, {a, (3, 7} subject to the 
normalization conditions. The total energy of the system is written as 

E B = 4tt / drr 2 £ B (r). (26) 
Jo 

We find the following result for the energy density which differs from GHGU Eq. [5.3]. The differences can be traced 
to different results for coherent-pair matrix elements. Writing the quark Dirac spinor as 

*tW=(.M*)xi m &, (27) 



v(r)a ■ f 



the energy density is given by 



Wr) = \{%) 2 + \{o\r) - v 2 ) 2 mlUa(r) 

+3[u(r)(& + fv{r)) v{r)f r + ga(r)(u 2 (r) - v 2 (r))} 

+ x){{f r ? + ^4> 2 {r)) + {N* - x)d>l(r) (28) 
— aS(a + b)u{r)v{r)4>{r) 

+\ 2 [x 2 + 2xN v + 81(a 2 a 2 c 2 + (f3 2 + 7 2 )d 2 ]0 4 (r) 
+\ 2 (N 7I + x)(a 2 (r)-v 2 )<f ) 2 (r), 

where N„ is the average pion number (GHGU Eq. [5.12]), 

N„ = 9[ a V + (/3 2 + 7 2 )c 2 ], (29) 

and where S takes the following values for nucleon or delta quantum numbers, 

S N = (5/3 + 4a/2 7 )/V3 (30) 
S A = (2V2/3 + 5 7 )/V3. (31) 

The function <p p {r) is obtained from <f>(r) by the double folding, 

POO 

4> p (r) = / Gj(r,r')<f>(r')r' 2 dr', (32) 
Jo 

with 

2 f°° 

oj(r,r') = - dkk 2 w(k)j 1 (kr)j 1 (kr'). (33) 
n Jo 

For fixed a, (3, and 7 the stationary functional variations are expressed by 

6[ / drr 2 (£ B (r) - 3e(u 2 (r) + v 2 {r)) - 2/s# p (r))] = 0, (34) 
Jo 

where the Lagrangian parameter, k, enforces the pion normalization condition, 

8tt / 4){r)(j) p {r)r 2 dr = 1, (35) 
Jo 

and the Lagrangian parameter e fixes the quark normalization, 



/•■DC 



4tt / drr 2 (u 2 (r) +v 2 (r)) = 1. (36) 
Jo 

Minimizing the Hamiltonian yields the four nonlinear coupled differential equations, 
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P = -(gcr + e)v(r) - §a<5(a + b)gtf>(r)u(r), 

% = - j!«(r) - (ga(r) - e)u{r) + §a<5(a + b)g(j>{r)v{r), 

& = <f* + 3 9(u 2 (r) - v 2 {r)) + 2\ 2 (N V + x)tf(r)a(r) + X 2 (a 2 (r) - «/>(r), 

+^(l + ^)(a 2 (r)-^(r) 



(37) 



N 



2-(a + b)g5u(r)v(r)-$-<f> p (r), 



with eigenvalues e and k . These consist of two quark equations for u and v where a(r) and <f>(r) appear as potentials, 
and two Klein-Gordon equations with u{r)v{r) and (u 2 (r) — v 2 {r)) as source terms. The boundary conditions are for 
r -> 0, 



and for r — > oo, 



i; = — = = -p =0, 



[r(ff/ w - £) 1/2 + + £) ( " 1/2) ] «(r) - r( 5 A + ef'Mr) = 0, 
(2 + 2m 7T r + m 2 r 2 )(f>(r) + r(l + m n r)(j>'{r) — 0, 
ra / (r) + (tr(r)-/ w )(l + m <r r)=0 ) 



(38) 



(39) 



which has one sign in the first equation different from GHGU Eq. [5.18]. The field equations are solved for fixed 
coherence parameter, x, and fixed Fock-space parameters, {a, (3, 7}. Then the expectation value of the energy is 
minimized with respect to {a, (3, 7} by diagonalization of the "energy matrix" , 



(40) 











a 




a 


H a /3 








a 


= E 


a 

















Each H entry of the matrix is related to a corresponding density, E(r), as follows: 

/■ 00 

H aP = 4ir r 2 E a0 {r) dr, 
Jo 

and analogously for the other entries. The E a p(r) functions for a nucleon arc 

E aa (r) = Eq(t) + 18a 2 <p 2 p (r) + 9a 2 X 2 {2x + 9c 2 )0 4 (r) 
+9a 2 X 2 [a 2 (r) - is 2 ](f> 2 (r) , 

Ep {r) = Eo(r) + 18c 2 2 (r) + 9A 2 (2xc 2 + 9d 2 )0 4 (r) 
+9c 2 A 2 [a 2 (r) -v 2 ]tf{r), 

5 



V3' 



E a p(r) = -2g(a + b)0(r)u(r)v{r) 



E ai {r) = -2g{a + b)(f>(r)u(r)v(r) • 4 

£ 77 (r) = E pi3 (r), 
E Pl {r) = 0. 

The E a p(r) functions for a delta are the same except for 

E a p{r) = -2g(a + b)<p{r)u(r)v(r) 
E al {r) = -2g(a + b)0(r)u(r)v(r) 



(41) 

(42) 
(43) 

(44) 

(45) 

(46) 
(47) 



2V2 

7r 
5 

71' 



which expresses the difference in the (5-terms (Eqs. [32-33]) appropriate to nucleon and delta respectively. 
In the above expressions, E (r) is given by 
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Eo(r) = - 

+U a + X 2 x 2 (/} 4 {r) + 3ga(r)[u 2 (r) - v 2 (r)] 

+3kt(r)(— + -w(rj) — I'M — — xm <p (r) 
dr r dr 

+\ 2 x[<j 2 (r)-v 2 ]cj) 2 (r). (48) 
Substituting for Uq and using the field equations, this can be rewritten as 

Mr) = \W(r)f + A 2 *W + 3ga(r)(u(r) 2 v{r) 2 ) + ^>(r) 2 f 2 ) 2 + ^(a(r) 2 f 2 ) 

-mlU{a{r) - h) + \ 2 x(a(r) 2 - f 2 )cf>(r) 2 - 2xml<\>{r)\ (49) 

which shows an explicit rapid decay as r — > oo. 

We solve this set of equations in the same iterative manner as GHGU. The iteration procedure is implemented as 
follows. For fixed values of x and a, (3, 7 the above differential equations with the corresponding boundary conditions 
are solved by using the same numerical package (COLSYS fLl|) as used in the original GHGU paper. Then the 
energy matrix is diagonalized and the minimum eigenvector chosen. These solutions are then mixed with the previous 
solution and repeated until self-consistency is achieved. The procedure is started with an initial guess. We used the 
so-called "chiral circle" meson field forms M, but the actual starting point does not matter, provided the iterations 
converge. 



da 
dr 



- ^ 



^[a 2 (r)-, 2 ] 2 - m lfMr) 



V. NUCLEON PROPERTIES 



In this section we review the several nucleon observables which will be calculated from the solutions arising from 
the procedure described in the previous section, noting differences with GHGU where they occur. From the electro- 
magnetic current operator, 

t = 1>{\ + ~7s)$ + e 3a0 $ a d^^, (50) 
one derives the charge and magnetic moment densities for the neutron and proton, 

P -M = a 2 (u 2 + v 2 ) + P 2 [\{u 2 + v 2 ) + U<t> p ] 
ine 3 3 

+7 2[4 (u 2 + t; 2 ) _2 #pL (51) 

+1 2 [-Uu 2 + v 2 ) + Ih p }, (52) 



47re 81 

x 



3" '3' 
^(54a 2 + 2/3 2 + 7 2 + 32V2/3 7 ) 

(9a 2 + a;)(4/3 2 +7 2 )0 2 , (53) 



729a 



2 



' (9a 2 + a ; )(4/3 2 -|-7 2 )0 2 , (54) 



729a 1 



which differ from GHGU Eq. [6.6] and [6.7] for the magnetic moment densities. The axial-vector to vector coupling 
ratio is given by (GHGU Eq. [6.9]) 

^ = 2(NJ Z = It 3 = ~\ [d 3 r: [-^ySW + ^d z arr - $ d z &] : \NJ Z = ±T 3 = \), (55) 
gy 2 2 J 2 2 2 

from which wc find, 



7 



5A 

9v 



An 



drr z 



25 



-a 



32V2 



da 



+ ^ + ^7 2 + ~irh] (u 2 - v 2 /3) + ^a(3(a + b)—</> 



27 



27 



3%/3 



dr 



(56) 



which differs with GHGU Eq. [6.10] in the first term. Finally, the 7rNN coupling constant can be calculated from the 
pion field or the pion source term. Using the pion field form, one has 



2M 



N 



aj3{a + b) drr 3 (j)(r), 



(57) 



which is in agreement with GHGU Eq. [6.15]. Using the pion source term, one can obtain an alternative form of the 
7rNN coupling constant, 



.9ttNN 

2M N 



4-7T 



, ., ,10 , 10 „, 50 , 64\/2 
drr A g(—a 2 + — /3 2 + — Y + 



50 , 
8l' 



(3^j)u(r)v{r) 



-4:TrX 2 a(3 / drr 
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81 A ' 81 ' ' 81 
JL(a + b)(o 2 {r) - f^(r) + -^=(a 2 b + 2b 2 a + ac 2 )^{r) 



(58) 



which differs from GHGU Eq.[6.14] by a factor of 3 in the last term. 

The cr-term was not calculated in GHGU, but is an important quantity in that it measures the degree of chiral 
symmetry breaking. For the linear sigma model considered here this quantity is given by 



<JirN = Airf^-ml / drr 2 (a(r) - 



(59) 



For a review of this quantity see Ref. p8[ 



VI. NUMERICAL RESULTS 



As mentioned earlier, the chiral quark model of the nucleon has two free parameters once the pion mass and pion 
decay constant are fixed at their experimental values. One can choose the free parameters to be the cr-mass and the 
meson-nucleon coupling constant, g. For comparison with GHGU and to illustrate the systematics of the model, for 
the results to follow, we fix these two parameters to m a — 700 MeV and <?=5.00. 

First consider the coherence parameter. Figure 1 shows the baryon (nucleon and delta) energy as the coherence 
parameter is varied. As was similarly shown in GHGU's Fig. [1], there is a clear minimum in the region of x=l with a 
corresponding N n = .43. Table I shows the various energy contributions to the baryon mass. While Table II gives the 
values of the derived parameters determined by the minimization procedure. The values found are in good agreement 
with that found in GHGU (self-consistent case) even though a different coherence dependence (which shows up in the 
d-term term in Eq. [29]) was used. This is expected however since the fourth-power pion self- interaction energy makes 
a relatively small contribution to the total energy. The nucleon-delta mass difference is found to be about 150 MeV. 
Figure 2 shows the quark wave functions and the meson fields for the x=l case, and Figure 3 shows the pion field 
and its derived kinetic energy density function, 4> p {r). All of these quantities are little changed from those reported 
by GHGU. 

Next consider the nucleon physical properties. Table II shows the results for several nucleon observables compared 
with the projected Hedgehog model of Birsc || and with experiment. For those nucleon quantities calculated by 
GHGU our results are nearly identical despite the differences in several quantities noted previously. 

In comparing the predictions of the coherent-pair model to those of the Hedgehog model, we find that the quark 
contributions to each observable are roughly similar; however the pion contributions are significantly smaller, a feature 
consistently seen throughout. This is somewhat surprising given the fact that both approaches are attempting to solve 
the nucleon problem from the same starting model albeit with very different methods and resulting different model 
parameters. Quantities such as the sigma-commutator which depend on the sigma field are quite similar to that of the 
Hedgehog model, both giving a sigma commutator of roughly 90 MeV. The small pion contributions to the calculated 
observables seems to be the coherent-pair approximation's principal phenomenological short-coming. 

Consider the nucleon charge radius squared. For the Hedgehog model the pion's contribution to the proton's charge 
radius squared is roughly 40% that of the quarks; while for the coherent-pair approximation the pions contribution is 
only 4% that of the quarks. The small pionic contribution is compensated for by a slightly larger quark contribution 
leaving the total proton charge radius squared very close to that calculated in the Hedgehog model, but still about 
20% that of the experimental value. 
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Likewise for the magnetic moments. The quark contribution to the magnetic moments is about the same for both 
models, but for the Hedgehog model the pions make a contribution to the proton magnetic moment which is 65% that 
of the quarks while for the coherent-pair approximation the pionic contribution is only 12% that of the quarks. The 
smaller pionic contribution to the magnetic moments results in magnetic moments roughly 60% that of the empirical 
values. 

For the case of the axial vector coupling constant we find again that the quark contributions in the two models are 
similar, but the mesonic contribution in the coherent-pair model is only half that of the Hedgehog model. In this case 
this actually helps in that the Hedgehog model predicts too large a value for this quantity. 

Finally we consider the 7rNN coupling constant calculated in two ways first using the pion field itself Eq.[59] and 
secondly using the pion source term Eq. [60] . In the Hedgehog model Birse found that both methods gave roughly the 
same value, but in the coherent-pair approximation the two methods give very different results. This large difference 
in the two, supposedly equivalent, ways of calculating g^NN has been noted already in Fiolhais, et al. || who studied 
the generalized projected Hedgehog model and compared with other models including GHGU. In Fiolhais et al. the 
difference in g^NN calculated from the pion field and the pion source terms along with the value expected from the 
Goldberger-Treiman relation provide a virial measure of how well the self-consistency condition in the pion sector is 
met. The Fiolhais et al. relation for the fractional virial deviation is 

t NN - 9-kNN , n l 1,U8 9nNN ~ 9-kNN 

g? NN (1-08)5^ + <V 



where g^NN is the source-calculated value, Eq. [60] , <7^Af Af ^ field-calculated value, Eq. [59] , ff^jv^y is the average of 
these two, and g^NN is the value expected from the Goldberger-Treiman relation, 

GTR _M^9A , ^ 

9-kNN — f ■ \ GL ) 
in gv 

The factor of 1.08 accounts for the effects of explicit chiral symmetry breaking. Using the coherent-pair approximation 
results of GHGU, Table 3 of Fiolhais et al. reports a fractional virial violation of 173%. This should be compared 
with 51% for the mean-field Hedgehog model of Birse and Bannerjee and 6% for the generalized projected Hedgehog 
model of Fiolhais et al. Using the values in Table III, we find a fractional virial violation of 149% which is a little 
better than that previously reported, but still shows a clear and substantial problem with self-consistency in the pion 
sector. 

One must question why this should be so. That the two models should give such different pionic contributions is 
a bit puzzling as the starting model is the same and both employ a self-consistent mean-field type approach which 
minimizes the energy. In comparing the quark wavefunctions and meson field solutions in the two models one finds, 
not surprisingly, that the quark wavefunctions and sigma field are nearly identical in the two models, but, though 
the pion field function in the coherent-pair approximation is smaller by a factor of about two from that found in the 
Hedgehog model, this cannot be the source of the discrepancy since the pion field in the coherent-pair approximation 
is normalized but the magnitude of the pion field (treated as a classical field) in the Hedgehog model is determined 
by the source terms. The magnitude of the pionic contribution to any quantity in the coherent-pair approximation 
is determined by the Fock space coefficients, 7} and by the various coherent functions {a,b,c,d}. With the 

minimized solution around x ~ 1.0 the mean number of pions in the coherent pair approximation case is N„ = 0.43 
which apparently yields the small pionic contributions to the various baryon observables. A comparison of the 
magnitude of the pion field resulting from the Hedgehog model with the normalized field arising in the coherent-pair 
model would imply an effective number of pions about 3 times greater in the Hedgehog model. A preliminary search 
of the parameter space to see if some other parameter set may correct this deficiency was not successful. We have 
examined the model up to a coherence parameter of x=2.25. Larger coherence parameters do indeed give rise to 
greater pionic contributions to the various nucleon observables, but at the expense of greater energy as well. While 
the issue is being investigated further, we suspect the problem with self-consistency as revealed in the large fractional 
virial violation to be the principal difficulty. 

Another useful comparison is shown in Figs. ^ and ^| where the energy densities for various terms are shown for 
both the Hedgehog model and the coherent-pair approximation. Again there is a systematically smaller contribution 
to each term in the coherent-pair model than in the Hedgehog model. Yet the total energy of both solutions are 
roughly the same (1070 MeV and 1120 MeV respectively). 

VII. SUMMARY AND CONCLUSIONS 

In this work we re-examined a linear sigma model of the nucleon using quarks and sigma and pion mesons as the 
fundamental degrees of freedom. To solve this model we have employed the coherent-pair approximation following 







the work of Goeke, et al. correcting several errors. We solved the model using a Fock space ansatz treating the sigma 
field as a classical field and treating the pions as quantum fields using the coherent-pair approximation of Bolsterli 
. We neglect vacuum effects and center-of-mass corrections. Despite the several corrections to the work of Goeke, 
et al., the numerical solutions we find are very close to those presented in their original paper as are all the calculated 
nucleon observables. We find that the calculated nucleon observables are reasonably close to experiment, but, in the 
case of electromagnetic quantities such as charge radii or magnetic moments, the pionic contributions seem too small 
when compared to that of other chiral nucleon models such as the Hedgehog model of Birse and Bannerjee The 
origin of this difference is not fully understood but probably arises from the lack of self-consistency in pionic sector 
using this approach as was noted previously by Fiolhais, et al.. Therefore, at this stage we must concur with the 
rather disappointing conclusion of Goeke, et al., namely, that a better description of the pionic sector is required 
before a more satisfactory nucleon phenomenology can be expected in this model. 
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TABLE I. Energy contributions to coherent-pair nucleoli using g=5 and m CT =700 MeV. (All values are in MeV.) 
Self-consistent solution evaluated using a coherent parameter value x=l. The details of this solution are discussed in the 



text. 

Quantity Nucleon Delta 

Quark eigenenergy 150 219 

Quark kinetic energy 1124 975 

Sigma kinetic energy 304 268 

Pion kinetic energy 236 185 

Quark-meson interaction energy -675 -318 

Meson interaction energy 84 114 

Baryon Mass 1073 1224 

Nucleon-Delta mass difference 140 



TABLE II. Numerical values of various parameters for the coherent-pair nucleon with a coherence parameter of x=l and 
using g=5 and m CT =700 MeV after self-consistent minimization. (N/A means "not available".) 



Quantity This work Goeke, et al. [1] 

Mass (GeV) 1.073 1.08 

a .820 .82 

P .379 .38 

7 .429 .43 

Pion eigenvalue, k (GeV) -.173 N/A 

Quark eigenvalue, e (GeV) .150 .14 

N„ .43 N/A 
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TABLE III. Observables for the coherent-pair nucleon with a coherence parameter of x=l and using m CT =700 MeV and g 
= 5. Magnetic moments are in nuclear magnetons. Charge radius is in fm. For comparison the results from the projected 
Hedgehog model of Birse [5] is also presented. 







Coherent Pair 






Hedgehog |J 






Quantity 


Quark 


Meson 


Total 


Quark 


Meson 


Total 


Exp. 


(r^) ch-proton (fni^) 


0.533 


.023 


0.556 


0.39 


0.16 


0.55 


0.70 


y } ch — neutron (flTl ) 


0.019 


-.023 


-0.004 


0.09 


-0.16 


-0.070 


-0.12 


(m) proton 


1.53 


.18 


1.71 


1.74 


1.13 


2.88 


2.79 


{/^) neutron 


-1.13 


-.18 


-1.31 


-1.16 


-1.13 


-2.29 


-1.91 


9A 

ff.iViv#Eq.[59]) 


1.07 


.39 


1.46 


1.11 


.75 


1.86 


1.25 






0.25 






0.93 


1.0 


^iVivf5f(Eq.[60]) 


1.11 


.24 


1.35 


1.16 


.379 


1.53 


1.0 


Sigma term (MeV) 






88.9 






94. 


~45 ±5 (l| 
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Nucleon and Delta Mass 

Coherence Parameter Dependence 



FIG. 1. Dependence of chiral nucleon and delta mass with respect to the coherence parameter, x, using g=5 and m CT =700 
MeV. 



Quark Wavefunctions and Meson Fields 
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FIG. 2. Quark wavefunctions and meson fields using coherence parameter of x=l with g—5 and m CT =700 MeV. 
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Pion Fields 
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FIG. 3. Pion field shown with kinetic energy pion field function, $ P (r), using coherence parameter of x=l with <?=5 and 
m CT =700 MeV. 



Coherent Pair and Hedgehog Models 
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FIG. 4. Comparison of kinetic energy densities for the Hedgehog model and the coherent-pair approximation using coherence 
parameter of x=l with g—5 and m a =700 MeV. 
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Coherent Pair and Hedgehog Models 

Quark-meson Energy Densities 
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FIG. 5. Comparison of quark-meson interaction energy densities for the Hedgehog model and the coherent-pair approximation 
using coherence parameter of x—1 with g=5 and m CT =700 MeV. 
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